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RADAR  IMAGING  FROM  A  DISTORTED  ARRAY 


INTRODUCTION 

A  microwave  array  may  be  distorted  for  many  reasons:  It  may  be  too  large 
to  be  surveyed  properly.  It  may  suffer  from  windloading.  Its  installation 
may  be  faulty.  Or  it.-?  distortion  may  be  electrical  rather  than  geometric: 
Medium  turbulence  can  cause  che  integral  of  the  dielectric  constant  over  a 
path  from  a  source  or  a  target  to  the  array  to  vary  randomly  with  position  in 
the  array.  Electromagnetic  coupling  from  the  antenna  elements  to  the  local 
environment  also  may  vary  randomly  with  position  in  the  array,  causing  random 
errors  in  the  driving  point  impedances  of  the  antenna  elements.  The  tolerance 
on  the  random  variation  is  about  30°  rms  or  somewhat  less  than  one-tenth 
wavelength  [1], 

That  a  microwave  array  designed  for  high  angular  resolution  imaging  is 
likely  to  suffer  from  one  or  more  of  these  difficulties  is  evident  from  an 
examination  of  (1) 

s(u)  ■  ^*i(x)  (/s(u)e^kxujdu]e  (1) 

L  j 

which  is  the  simplest  form  of  the  integral  'equation  relating  a  source  function 
or  scene  s  and  its  image  or  estimate  s  wherj  the  radiation  field  due  to  the 
source  is  measured  by  a  line  aperture  of  extent  L  having  aperture  weighting  i. 
The  source  is  assumed  to  be  in  the  far  field  of  the  aperture,  which  permits 
its  description  in  terms  of  the  one-dimensional  reduced  angular  variable 
u  ■  sinQ,  9  being  measured  from  the  normal  to  the  aperture.  The  aperture  also 
is  assumed  to  be  one-dimensionai.  The  exponential  kernels  are  Fourier  kernels 
and  both  integrals  are  Fourier  integrals.  The  inner  integral  is  the  radiation 
field  S  at  the  aperture.  (1)  can  be  written  in  terms  of  S: 

s(u)  -  J'  i(x)S(x)e  ^xudx 

L  (2) 

U(x)S(x)> 

where.?  {•}  means  Fourier  transform.  By  the  multiplication-convolution 
property  of  the  Fourier  transform 


s(u')  ■  J"  {i(x)}*3"{S(x) } 
■  I(u)*s(u) 

Manutcript  submitted  August  20,  1981. 


(3) 
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where  *  means  convolution.  I(u)  is  the  radiation  or  diffraction  pattern  of 
the  aperture.  Since  convolution  hroadens  or  spreads  a  function,  the  closer 
that  I(u)  approximates  a  5-function,  the  better  the  estimate  s(u)  is  of  s(u). 
Since  the  width  of  I(u)  is  always  the  order  of  X/L,  it  is  evident  that  the 
larger  the  aperture  size  L  the  more  closely  will  s  be  represented  by  s. 

The  same  conclusion  pertains  to  a  target  or  scene  in  the  near  field. 

-4  -5 

To  achieve  a  resolving  power  of  10  to  10  rad,  which  is  typical  of 
common  optical  instruments  such  as  cameras  and  small  telescopes,  a  microwave 
antenna  must  be  hundreds  of  meters  to  tens  of  kilometers  in  size.  Such 
antennas  are  too  large  to  be  constructed  as  single  structures  such  as  the 
parabolic  dish.  Instead  they  must  be  phased  arrays.  They  must  be  highly 
thinned  (mean  interelement  spacing  >>  X/2)  to  control  costs,  and  the  element 
distribution  must  be  aperiodic  to  eliminate  grating  lobes  [1],  It  is  most 
likely  that  the  large  size  will  preclude  accurate  knowledge  of  element  loca¬ 
tion;  hence,  the  array  properties  will  tend  to  degenerate  to  those  of  the 
random  array  and  the  system  design  must  accomodate  the  poor  sidelobe  perform¬ 
ance  expected  from  such  an  array  [2],  As  a  consequence,  array  distortion  is 
a  highly  likely  property  of  a  huge  array  designed  for  microwave  imaging. 
However,  its  effect  can  be  neutralized  to  a  considerable  extent  by  adaptively 
controlled  phase  corrections  within  the  system  based  upon  measurements  at  the 
array  of  the  radiation  field  from  a  point  source  [3,4].  Adaptive,  retrodirec- 
tive  beamforming  techniques  then  focus  the  distorted  array  upon  the  source 
[5,6]. 

The  paper  describes  an  algorithm  for  phase  synchronizing  or  self-cohering 
a  distorted  array  upon  such  a  source,  moving  the  focused  beam  in  range  and 
angle  to  a  target  area,  and  scanning  it  across  the  target  to  image  it.  An 
instrument  embodying  this  procedure  is  called  a  radio  camera.  Experimental 
evidence  of  the  validity  of  the  technique  is  given.  The  experiments  were 
conducted  with  an  X-band  (X  ■  3  cm)  radar,  using  a  distorted,  quasi-linear , 
random  antenna  array  27  m  long. 


THE  RADIO  CAMERA 


Figure  1  shows  a  badly  distorted  receiving  array  and  a  point-source  of 
radiation  in  its  near  field.  It  also  shows  a  pulsed  transmitter.  The  near¬ 
field  source  may  be  active,  such  as  a  beacon,  or  passive,  such  as  a  large 
cross-section  reflector  echoing  the  radiation  from  the  transmitter.  This 
source  or  reflector  is  called  the  adaptive  beamformer;  it  also  is  called  the 
phase  synchronizer.  The  arrav  ic  shown  measuring  the  phase  of  the  signal 
received  at  e.  ch  element  relative  to  che  phase  at  some  reference  element 
within  the  array.  In  the  absence  of  multipath  the  phase  $n  of  the  received 

signal  at  the  nth  element  is  (tor  /c  )  +  C  ,  where  to  is  the  radiation  frequency 

n  n  n 

r^  is  the  path  length,  c^  is  the  average  speed  of  propagation  over  the  path 

and  $n  is  the  phase  shift  through  the  receiving  element  and  associated 

circuits,  -i  mav  be  rewritten  as 
n 


D 


n 


w(r  +  5r  ) 
on _ tv_ 

c  +  dc 
o  n 


+  ’n 


(4) 
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Fig.  1  —  Distorted  array  self -coheres  on  a  phase  synchronizing  source 


to  explicitly  show  the  displacement  6r^  of  the  element  in  the  direction  toward 

the  adaptive  beamformer  from  its  correct  distance  r  ,  and  the  deviation  5c 

on  n  ' 

of  the  mean  propagation  speed  over  the  nth  path  from  the  average  speed  c 

in  the  medium.  Since  6c  <<  c  , 

n  o 


w(r  +  6r  ) (1  -  6c  /c  ) 

*  *  — SSL - 2 _ B— g.  4  e 

n  c  n 

o 


kr  5  c 

=  kr  +  k6r - 22 — 2  +  5  k  -  — 

on  n  c  n’  c 


The  phase  differences  from  element  to  element  are  seen  to  have  four 

components.  The  primary  component  kr  is  due  to  source  location  relative  to 

on 

the  undistorted  array  and  is  the  same  as  in  any  phased  array.  The  residuals 
are  due  to  the  distorted  geometry  of  the  array,  to  variations  in  propagation 
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conditions  from  the  source  to  different  parts  of  the  array  due  to  its  large 
size,  and  to  variations  in  electromagnetic  coupling  between  the  elements  and 
their  surroundings.  By  phase-shifting  each  receiver  channel  by  the  negative 
of  these  phase  differences,  all  signals  from  the  adaptive  beamformer  become 
cophased.  The  output,  following  summation,  is  that  of  an  array  focused  upon 
the  beamforming  source.  The  spot  size  of  the  focal  zone  in  the  near  field, 
calculated  from  diffraction  theory,  is  shown  in  the  figure.  It  is  nominally 
XR/L  in  the  cross-range  dimension  and  7X(R/L)2  in  range.  X/L  is  the  beamwidth 
of  a  focused  diffraction-limited  aperture  of  length  L;  it  also  is  the  beam- 
width  of  the  distorted  array  following  the  adaptive  beamforming  procedure. 
7X(R/L)2  is  the  approximate  value  of  the  3  dB  range  beamwidth  due  to  the 
limited  depth  of  field  of  an  aperture  focused  on  a  near-field  target  [1]. 

XR/L  also  is  the  far-field  cross-range  dimension. 

Also  shown  schematically  in  the  figure  is  a  second  set  of  phase  shifters 
located  prior  to  the  summer.  Unlike  the  first  phase  shifting  operation,  the 
second  operation  is  open  loop  and  nonaaaptive.  After  the  array  is  phase  syn¬ 
chronized  the  focused  beam  is  scanned  in  range  and  angle  by  open-loop  correc¬ 
tions  calculated  from  the  geometry.  The  calculations  are  made  exactly  as  in 
a  conventional  phased  array.  During  near-field  scanning  the  spot  size  and 
shape  remain  as  indicated  in  the  figure.  In  the  very  near  field  the  resolu¬ 
tion  of  the  two-dimensional  image  that  results  from  the  scanning  operation 
is  about  equal  to  the  spot  size.  When  the  target  is  in  the  far  field  of  the 
large  array  the  depth  of  field  becomes  infinite  and  the  resolution  in  range 
is  determined  by  the  pulse  duration  of  the  transmitter.  At  intermediate 
distances  the  range  resolution  is  the  smaller  of  the  depth  of  field  and  the 
radar  pulse  length. 

Both  phase  shift  operations  can  be  analog  or  digital.  The  phaselock  loop 
is  the  natural  analog  circuit  for  adaptive  beamforming.  The  experimental 
equipments  described  later  use  digital  phase  shifting.  Figure  2  shows  the 
procedure  more  explicitly.  A  transmitted  pulse  illuminates  both  a  target 
area  to  be  imaged  and  a  passive  phase  synchronizing  target,  which  is  sketched 
as  a  corner  reflector.  The  echo  trace  consists  of  their  echoes  plus  clutter. 
The  signals  received  at  the  several  antenna  elements  are  sampled  in  range  and 
stored  in  the  format  shown.  The  range  trace  delivered  by  each  antenna  element 
is  stored  as  a  row  of  complex  numbers.  Successive  rows  correspond  to  succes¬ 
sive  elements  in  the  array.  The  position  of  a  sample  in  a  row  is  proportional 
to  range  and  designates  the  range  bin.  The  position  of  a  sample  in  a  column 
designates  the  array  element  number  and  is  monotonic  with  but  not  necessarily 
proportional  to  element  position  since  the  array  is  distorted. 

For  simplicity  it  is  assumed  that  ! r  -r  I  <<  A R,  all  n,  m,  where  AR  is 

n  m 

the  length  of  a  range  cell.  Thus,  all  echoes  from  a  common  target  will  appear 
in  a  single  range  bin  in  the  format  shown  in  Figure  2.  Actually,  target  echoes 
can  appear  in  different  range  bins  for  two  reasons.  The  first  is  due  to  the 
large  size  of  the  array,  for  when  Lsin  0  >  AR,  the  differential  arrival  time 
from  a  source  at  angle  9  exceeds  the  pulse  length.  This  effect  is  easily 
calculated  and  therefore  may  be  corrected  in  the  signal  processor.  The  aecond 
cause  is  array  distortion  in  the  direction  toward  the  source.  Most  arrays 
.designed  to  be  linear  or  planar  systems  will  not  suffer  distortion  so  severe 
as  to  warrant  range-bin  correction,  yet  will  still  require  adaptive  phase  cor¬ 
rection  if  the  unknown  a  priori  geometric  distortion  exceeds  a  small  fraction 


of  a  wavelength.  Large  irregular  arrays  having  deviations  from  planarity 
comparable  to  or  exceeding  AR  also  will  require  range-bin  correction. 

Assuming  that  range-bin  correction,  if  needed,  has  been  accomplished, 
three  distinct  additional  operations  upon  the  stored  data  are  required  for 
imaging.  These  constitute  the  algorithm  and  are  given  in  separate  sections 
below.  Two  are  uhe  beamforming  and  scanning  operations  discussed  earlier. 
The  third,  which  precedes  these  two,  is  the  search  procedure  for  the  phase 
synchronizing  source. 
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REFERENCE  SIGNAL 

The  first  operation  of  the  signal  processor  following  data  sampling  and 
storage  is  a  search  in  range  for  a  good  phase  synchronizing  source,  i.e.,  for 
that  echo  sequence  across  the  array  that  most  closely  approximates  the  expected 
field  from  a  point  source.  That  range  is  designated  R  and  is  the  reference 
range,  i.e.,  the  range  from  which  the  phase  synchronizing  reference  signal  is 
obtained.  The  radiation  field  from  a  point  source  in  free  space  would  be  near¬ 
ly  constant  in  amplitude  across  the  array,  while  the  phase  differences  would 
disclose  the  differential  distances  (modulo  wavelength)  to  the  source  as  well 
as  the  differential  phase  shifts  through  the  receiving  elements.  A  simple  test 
to  find  that  target  whose  radiation  field  most  closely  approximates  that  of  a 
point  source  is  to  measure  the  normalized  echo  amplitude  variance  at  each 
range.  R0  is  that  range  for  which  che  test  value  is  minimum. 

Figure  3  shows  typical  echo  amplitude  sequences  measured  across  a  large 
array.  One  is  from  a  corner  reflector  and  the  other  from  a  pick-up  truck  with 
camper  top,  which  is  a  complicated  target.  The  array  was  a  100  sample  point, 
distorted  aperture  27  m  in  length  operating  at  3  cm  wavelength  [3].  The  corner 
reflector  is  the  ideal  adaptive  beamformer,  yet  the  echo  amplitude  pattern  is 
not  constant  with  element  position.  There  are  four  contributing  factors  to 
the  general  variation  of  echo  amplitude  with  element  position.  The  first  is 
due  to  the  size  of  the  reflector  or  source.  The  second  is  due  to  multipath. 

The  third  is  due  to  clutter.  The  fourth  is  due  to  the  radiation  pattern  of 
the  antenna  element  used  in  the  array. 


Fig.  3  —  Echo  strength  from  two  range  bins  vs  element  position  in  array. 
The  comer  reflector  approximates  a  point  source. 
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A  reflector  of  size  T  reradiates  a  lohular  pattern  in  which  the  nominal 
lobe  spacing  is  X/T  rad  (Figure  4,  left).  The  nominal  lobe  spacing  at  the 
array,  which  is  at  a  distance  R,  is  XR/T,  and  the  echo  sequence  across  the 
array  would  be  expected  to  have  a  spatial  period  of  this  value.  If  XR/T  <  L 
a  low  spatial  frequency  or  long  period  amplitude  modulation  would  be  e-'ident 
in  the  measured  data.  Such  a  long  period  spatial  modulation  appears  in  the 
corner  reflector  echo  sequence  in  Figure  3  (although  the  cause  is  not  due  to 
the  size  of  the  reflector  but  to  the  radiation  pattern  of  the  antenna  element, 
which  is  discussed  later). 

The  high  spatial  frequency  or  short  period  fluctuations  seen  in  the  corner 
reflector  echo  sequence  in  Figure  3  are  due  to  clutter  (Figure  4,  middle).  The 
transmitter  beainwidth  is  broad  compared  to  the  high-resolution  receiver  beam- 
width;  the  beamwidth  A9  =  X/a,  where  a  is  the  size  of  the  radiating  antenna. 

The  clutter  paten  illuminated  by  the  transmitter  at  distance  R  is  RA0  =  RX/a. 

The  lobes  of  the  backscattered  radiation  have  a  nominal  spacing  X/(RXa)  *  a/R, 
and  the  nominal  cross  section  at  the  array  is  equal  to  R(a/R)  =  a.  This  last 
expression  means  that  the  correlation  distance  of  the  spatial  amplitude  modu¬ 
lation  due  to  the  clutter  equals  the  size  of  the  transmitting  antenna,  which, 
in  this  experiment,  was  less  than  the  average  interelement  distance.  Since  the 
transmitting  antenna  is  small  compared  to  the  large  receiving  array,  the  clutter 
modulation  always  will  appear,  as  in  Figure  3,  as  a  high  spatial  frequency  or 
short  period  fluctuation. 

The  effect  of  multipath  is  similar.  Let  $  be  the  angular  separation  at  the 
array  between  the  direct  arrival  and  the  multipath  signal.  The  complex  ampli¬ 
tude  of  the  received  signal  along  the  array  (x  coordinate)  has  the  form 
1.+  mexp(jkxS)  where  a  is  the  reflection  coefficient  of  the  multipath  scatterer 
and  k  =  2ir/X  is  the  wavenumber.  The  period  of  the  amplitude  modulation  of  the 
sum  of  the  two  signals  is  X/g.  Rarely  is  8  >  A6 ;  therefore,  multipath  in 
general  will  introduce  a  lower  spatial  frequency  into  the  echo  sequence  across 
the  array. 

The  fourth  contributing  factor  to  the  nonconstancy  of  the  reference  echo 
sequence  is  the  radiation  pattern  of  the  receiving  element  used  in  the  array. 
Figure  5  pictures  four  array  elements  each  of  length  d  in  an  array  of  length  L. 
The  beamwidth  of  an  element  pattern  is  X/d  and  its  cross  section  at  the  beam¬ 
forming  target  is  XR/d.  The  figure  illustrates  that  when  L  >  XR/d  the  element 
gain  to  the  target,  varies  with  element  position  in  the  array. 

The  echo  variance  calculation  should  be  based  upon  the  first,  second  and 
third  factors,  which  measure  the  quality  of  the  radiation  field  for  adaptive 
beamforming,  but  not  the  fourth  factor,  in  which  the  measuring  instrument 
induces  a  variation.  The  effect  of  the  latter  is  minimized  by  dividing  , 
the  echo  amplitude  from  the  ith  range  bin  received  by  the  nth  element, 
by  the  estimated  element  pattern  gain  f  from  the  nth  element  to  the  target. 

(A  single  subscript  suffices  for  the  element  pattern  gain  unless  the  target 
is  in  near  field  of  the  element,  which  is  highly  unlikely.) 
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Fig.  4  —  Target  size  and  clutter  contribute  to  echo  amplitude  variation  across  array 


The  corrected  amplitude  is 


*“  _  1  2 
A,  =  V,  /f  .  Its  mean  A,  *  -  I  A,  and  mean  square  At  = 
in  in  n  i  N  ,  in  ^  i 

n=l 

are  calculated,  where  N  is  the  number  cf  antenna  elements. 


1  N  2 
=.  1 

N  .  in 
n=l 

Its  variance 


2  1  _  2  2  _  2 

o,  =  -  I  (A.  -  A.)  =  A.  -  A.  .  o.  is  normalized  to  A 

i  N  ,  in  i  i  i  x 

n-“l 

effect  of  target  strength  and  range  from  the  calculations. 


2 

1 


to  remove  the 
The  normalized 


—  2 

variance  is  1  -  A^  /A^  ,  which  is  a  minimum  at  that  range  element  for  which 
— ^  2 

Ai  ^Ai  larSest  value.  The  rule  for  finding  Rq,  then,  is  to  calculate 

Ai  f°r  ran8e  elements  and  search  for  the  largest  value.  The  signal 

sequence  from  that  range  bin  becomes  the  reference  signal. 
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ADAPTIVE  BEAMFORMING 


Table  1  shows  all  the  steps  in  the  procedure.  Step  1  initiates  the  process. 
V^expCj’^^)  is  the  complex  envelope  of  the  echo  from  the  ith  range  bin  re¬ 
ceived  by  the  nth  element.  Steps  2  and  3  are  the  search  for  R0  discussed  above. 
Step  k  is  the  adaptive  beamforming  step.  The  measured  phases  ^  can  be  broken 

into  the  sum  of  two  terms  iJj!  +  4»! '  .  The  first  term  contains  the  conventional 

m  in 

target  and  array  geometry  and  is  all  that  would  be  expected  in  the  absence  of 
array  distortion,  medium  turbulence,  multipath  and  scattering.  The  second  term 
represents  the  errors  which  the  phase  synchronization  process  must  overcome. 

In  the  absence,  of  such  errors  the  image  or  the  target  or  clutter  from  the  ith 
range  element  would  be  obtained  from  the  integral  of  the  phase-weighted  signals 
received  acres.?,  the  array  from  the  ith  range  bin.  As  in  any  phased  array  the 
phase  weighting  is  the  conjugate  of  the  kernel  of  the  diffraction  integral. 

For  ease  of  discussion  the  diffraction  integral  can  be  approximated  by  the 
Fresnel  integral,  which  reduces  to  the  Fourier  integral  when  the  target  is  in 
the  far  field.  That  integral  in  turn  is  best  represented  by  a  sum  as  in  (6) 
because  the  array  is  discrete.  The  sum 


s.(u) 


N 

Z 

n-1 


1^.  -]k(x  u 

.  J  m  J  n 
A,  e  e 
ir. 


V2V 


(6) 


is  the  image  s^  of  the  scattering  sources  in  the  ith  range  bin  provided  that 
i|)M  =0.  As  in  (1),  u  =  sin  8  and  0  *  scan  angle  from  the  normal  to  the  array. 
The  discrete  variables  xn  and  R^  are,  respectively,  the  x  coordinate  in  the 
array  (transverse  to  the  array  normal)  of  the  nth  element  and  the  distance  to 
the  ith  range  bin.  Although  s^Cu)  is  written  as  a  one-dimensional  image 
(in  the  reduced  angular  variable  u)  of  the  echoes  from  the  ith  bin,  it  is,  in 
reality,  a  two-dimensional  image  in  u  and  R.  Consequently,  two-dimensional 
imagery  is  available  for  display  of  near-field  targets. 


The  function  of  the  adaptive  processor  is  to  compensate  for  so  that 
the  operation  described  by  (6)  may  be  accomplished.  The  earlier  steps  in  the 
process  searched  the  data  to  find  the  range  bin  in  which  the  echo  amplitude 
across  the  array  had  the  smallest  variation.  That  range  bin  became  the  refer¬ 
ence  range  and  the  signal  from  that  range  bin  became  the  synchronizing  signal. 
Its  complex  envelope  following  normalization  to  f  is  Aonexp(j’4>on)  where  the 
amplitudes  are  all  nearly  the  same.  The  phases  of  the  echoes  are  raiidom, 
however,  due  to  the  perturbed  geometry  of  the  array  or  the  spatial  variations 
of  the  refractive  index  or  the  impedance  variations  from  element  to  element  or 
any  combination  of  these  factors.  If  the  reference  source  were  an  ideal  source 
the  phase  differences  would  be  due  entirely  to  them.  In  addition  the  phases 
are  perturbed  by  multipath  and  scattering  and  by  the  finite  lobe  width  of  the 
reradiation  from  the  synchronizing  source.  However,  because  it  is  generally 
impossible  for  the  system  designer  to  obtain  a  priori  information  about  these 
latter  conditions,  the  signal  processor  must  necessarily  ignore  them. 

The  fourth  step  in  the  process  is  to  compensate  for  the  phase  variations, 
which  are  assumed  to  be  due  to  the  first  set  of  factors.  Correction  is  accom¬ 
plished  by  phase  rotating  the  complex  envelopes  of  the  signals  from  Rq  received 
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by  the  different  array  elements.  The  proper  phase  shift  for  the  nth  element 
is  the  negative  of  the  phase  difference  A<|>on  *  -  <{i  .  The  complex  signal 

envelope  at  the  nth  element  becomes  Aonexp(j^on)exp[-j  (<J^n  -  ^00)]  3  Aexp(jC'  ). 

This  correction,  or  phase  conjugation,  is  exactly  what  a  phased  array  or  lens 
would  do  when  focusing  upon  the  reference  source.  The  output  signal  from  the 
array  when  it  is  so  focused  is  the  sum  of  these  phase-corrected  echoes  from  R  . 


Table  1  —  Steps  in  Radio  Camera  Imaging 
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SCANNING 


The  fifth  step  (performed  simultaneously  with  Step  4)  is  to  phase  rotate 
all  the  signal  samples  from  each  antenna  element  by  i{ion  -  d»00 . 

The  samples  of  the  complex  envelope  from  the  ith  range  bin  now  become 
Alne*p[.H*1„  -  *on  +  *„„>]•  The  sum  of  such  a  set  of  complex  samples  represents 

the  output  of  a  misfocused  phased  array  or  lens  since  iK  4  i>  except  at  the 

range  R  of  the  reference  reflector.  n  on 

o 

The  sixth  step  is  to  focus  the  array  at  all  ranges  simultaneously.  This 
task  is  accomplished  for  an  arbitrary  range  by  refocusing  the  array  from 
the  reference  range  Rq  to  R^.  The  phase  correction  is  approximately  quadratic, 
as  given  by  (6).  Assuming  that  the  earlier  steps  were  performed  properly,  the 
self-cohering  process  forced  the  quadratic  component  of  the  phase  shift  of  the 

2  2 

signal  in  the  nth  channel  to  become  -kx  /2R.  +  kx  /2R  .  To  focus  the  array  to 

n  i  n  o 

range  this  term  must  be  set  to  zero,  which  requires  a  further  phase  addition 

cf  (kx^/2) (1/R^  -  1/Rq) .  This  step  requires  a  knowledge  of  the  range  RQ  of  the 

reference  reflector.  Fortunately,  the  value  of  Ro  is  available  in  the  system 
for  ic  is  measured,  as  in  conventional  radar,  by  the  round-trip  travel  time 
of  the  pulse  to  the  phase  synchronizing  source  and  it  is  read  directly  into 
the  signal  processor  from  the  radar  receiver.  The  accuracy  of  measurement  is 
determined  by  the  range  resolution  of  the  system,  which  is  the  order  of  the 
reciprocal  of  the  signal  bandwidth  (in  distance  units) ,  or  the  near-field  range 
beamwidt.h,  whichever  is  smaller. 

Step  7  imparts  a  linear  phase  rotation  to  the  range-focused  complex  envelope 
(designated  B^  in  Table  1)  for  each  scan  angle  u.  The  phase  shift  is  -kxnu. 

The  last  step  forms  the  sum  of  the  linearly  phase-weighted,  range-focused 
samples  to  obtain  the  image  s^(u)  for  the  ith  range  bin: 


^(u) 


N 

Z 

n-1 


B.  e 
in 


-jkx  u 
J  n 


(7) 


ELEMENT  POSITION  TOLERANCE 


Steps  3,  4,  and  5  require  no  knowledge  whatsoever  of  element  position. 

The  phase  synchronization  process  is  purely  retrodirective.  Steps  6  (refocus¬ 
ing  in  range)  and  7  (scanning  in  angle)  do  require  coordinate  information. 

The  tolerances  on  element  position  error  have  been  worked  out  [7],  [8], 

[1,  Chapter  13].  The  most  stringent  tolerance  is  invoked  by  Step  7.  Briefly, 
the  theory  consists  of  the  following  points.  First,  the  loss  in  main  lobe 
gain,  in  dB,  due  to  all  the  random  phase  errors  across  the  array  is 

i  2 

LG  -  4.3c“  where  n  i3  the  variance  of  the  phase  errors  in  square  radians. 


Second,  the  phase  variance  due  to  random  position  errors  is  approximately 
2  2  2  2  2 

=  k  where  a*  is  the  position  error  variance  in  the  array  in  the 

direction  perpendicular  to  the  beamforming  direction  and  6  is  the  scan  angle 
measured  from  the  direction  of  beamformina. 
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Combining  these  expressions,  and  allowing  one  dB  loss  in  array  gain  as  an 
acceptable  gain-loss  tolerance,  the  rms  element  tolerance  becomes 


The  factor  X/4it  is  the  conventional  tolerance  in  phased  array,  mirror  or 
lens  design.  9max(rad)  is  the  maximum  scan  angle  from  the  direction  of  the 
reference  source  and  is  half  the  field  of  view.  Since  an  individual  radar 
target  always  subtends  a  very  small  angle  at  the  radar,  0max  <<  1  and  the 
allowed  element  position  error  is  exceedingly  large.  For  example,  if  the  sys¬ 
tem  were  designed  to  image  a  target  as  large  (in  angle)  as  the  moon,  @max 
would  be  approximately  10"2  rad;  in  this  case  the  position  accuracy  tolerance 
increases  by  two  orders  of  magnitude  to  about  10X .  It  is  precisely  this 
extraordinary  liberty  in  position  tolerance,  following  adaptive  phase 
synchronization,  that  permits  radio  camera  imaging  with  a  nonrigid  or  ill- 
surveyed  array. 


EXPERIMENTS 

Figure  6  shows  the  result  of  an  X-band  radio  camera  imaging  experiment 
using  a  low  power  (5w)  pulsed  radar  transmitter  and  a  corner  reflector  for  the 
phase  synchronizer.  The  size  of  the  reflector  was  0.46m  and  its  radar  cross- 
section  was  calculated  to  be  56m2.  The  estimated  clutter  cross-section  was 
2m.  The  receiving  array  was  the  one  use  in  the  experiment  for  Figure  3;  it 
was  27m  long  and  1.2m  wide.  The  100  element  positions  were  randomly  located 
within  it.  Uniform  probability  density  functions  of  element  position  were 
chosen  for  both  the  length  and  width  dimensions  of  the  array.  The  experiment 
was  conducted  on  a  time-shared  basis  in  which  a  single  radar  receiver  was 
successively  moved  from  position  to  position,  delivering  a  radar  echo  trace  to 
a  microprocessor  from  each  radar  position,  after  which  the  operations  described 
earlier  were  performed.  The  image  is  one-dimensional,  deflection  modulated. 

The  target  consisted  of  two  additional  reflectors,  each  0.61m,  drawn  to  scale 
below  the  image.  Target  distance  from  the  array  was  240m.  The  reference 
reflector  was  41m  closet.  Drawn  also  in  Figure  6  ^s  the  calculated  response 
of  the  array  in  free  space  had  it  self-organized  perfectly.  It  is  evident 
that  the  experiment  was  exceedingly  successful. 

The  theoretical  beamwidth  for  this  experiment,  based  upon  diffraction 
theory,  was  0.88X/Lcos9.  The  coefficient  corresponds  to  the  particular 
probability  density  function  of  element  location  used  [1].  0,  the  target 

angle  from  the  array  normal,  was  25®.  This  expression  evaluates  to  l.lmrad, 
which  is  indistinguishable  from  the  measured  beamwidth. 

Figure  7  is  the  image  of  the  same  reflectors  when  the  equipment  was  oper¬ 
ated  in  the  synthetic  aperture  mode  [9]:  the  low  power  transmitter  and  receiver 
both  were  moved  from  position  to  position  for  each  radar  transmission  and 
reception.  (The  experimental  set  and  procedure  were  illustrated  and  described 
in  [3].)  Because  of  the  doubling  of  the  wavenumber  due  to  the  synthetic  aper¬ 
ture  operation,  the  beam  cross-section  in  Figure  7  is  halved.  Again,  the 
comparison  with  the  calculated  free  space  response  is  excellent. 

In  both  experiments  the  sidelobe  properties  conform  to  the  theory  of  the 
random  array  [1,2],  which  predicts  that  the  contribution  to  the  average 
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Fig.  6  —  One-dimensional  image  of  two  comer-reflector  targets  using  another 
comer  reflector  for  adaptive  beamforming.  Theoretical  response  in  free  space 
shown  for  comparison. 


14 


30  f- 


SIMULATED  TARGETS 


'3 


1 


*  TARGETS 


U—  0.61  m  (SUBTENDS  2.54  mr 
1  AT  TARGET  RANGE) 


Fig.  7  —  Radio  camera  image  of  same  two  targets  in  synthetic  aperture  mode 
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sidelobe  power  level  is  N~^  times  the  main  lobe  power  response  of  a  target. 
Thus,  the  theoretical  average  sidelobe  level  for  the  two  nearly  equal  strength 
targets  in  Figures  6  and  7  is  -17  dB;  the  measured  levels  are  within  a  few 
tenths  of  a  decibel  of  this  value. 


ISOLATED  TARGET 

The  situation  is  somewhat  different  for  an  isolated  target  such  as  an 
aircraft  or  a  ship  in  which  the  reference  source  is  on  the  target.  One  differ¬ 
ence  is  that  no  range  refocusing  (Step  6  of  Table  1)  is  required  since  the  size 
of  the  aircraft  usually  will  be  small  compared  to  the  depth  of  field  of  the 
array.  A  more  significant  difference  lies  in  the  fact  that  target  echoes  are 
used  both  for  adaptive  phase  synchronization  as  well  as  for  imaging.  The 
requirements  upon  the  echo  characteristics  are  opposed  for  these  two  processes. 
Phase  synchronization  requires  a  dominant  point  source,  whereas  the  objective 
of  picture  taking  is  to  reproduce  the  angular  backscatter  profile  of  a  compli¬ 
cated  target  so  that  it  may  be  recognized. 

The  conflict  resolves  itself  when  the  target  is  moving  relative  to  the 
observer.  Figure  8  shows  an  aircraft  at  distance  R  wich  its  velocity  vector  V 
making  an  angle  a  with  the  direction  to  the  array.  The  aspect  angle  a  changes 
with  time  at  the  rate  da/dt  «  V  sina/R.  Because  of  this  turning  motion,  the 
radar  echo  changes  with  time,  sometimes  having  the  desirable  properties  for 
phase  synchronization  while  at  other  times  the  return  signal  is  good  for  imag¬ 
ing  purposes.  From  time  to  time  a  strong  highlight  appears  due,  say,  to  a 
broadside  specular  return  from  the  fuselage,  or  a  two-plane  corner  reflector 
formed  by  wing  and  fuselage  or  in  the  tail  assembly.  The  highlights  are 
severely  aspect-angle  dependent,  as  is  demonstrated  in  Figure  9.  There  a  small 
flat-late  reflector  of  length  T  is  shown  momentarily  oriented  broadside  to  the 
direction  to  the  array.  As  the  line-of-sight  angle  6  changes  due  to  platform 
motion  the  reradiation  pattern  moves  through  twice  that  angle.  Consequently, 
the  time  that  the  main  lobe  dwells  upon  the  array  is  determined  by  the  turning 
rate  da/dt  and  the  reradiation  beamwidth  («X/T).  The  duration  or  correlation 
time  can  be  estimated  as  the  time  required  for  the  differential  phase  of  the 
echoes  from  two  scatterers  separated  a  distance  equal  to  the  target  size  T  to 
change  about  ir/2.  The  right  of  Figure  9  shows  the  geometry.  The  change  in  the 
differential  distance  SR  from  two  such  scatterers  is  approximately  Tda.  The 
rate  of  change  Is  d(5R)/dt  3  Tda/dt  3  TVsina/R  and  the  phase  difference  between 
their  echoes  changes  at  a  rate  i(<5$)/dt  *  2kd(6R)dt  3  unTVsina/XR.  Hence,  the 
correlation  time  is  the  order  of 


T  3  XR/8TVsina  (8) 

c 

Consider  a  Mach  1  aircraft  with  a  reflecting  region  one  meter  in  size 
oriented  normal  to  the  radar  line-of-sight.  When  observed  with  3  cm  radiation 
from  a  distance  of  20  km  and  at  a  309  angle  the  correlation  time  is  about  1/2  s. 
During  this  time  the  large  array  on  the  ground  may  be  synchronized,  forming  a 
narrow  beam  focused  on  the  specular  reflector.  Rather  than  determining  at  what 
range  to  phase  conjugate  or  focus,  as  in  the  earlier  discussion  (Step  3  of 
Table  1),  the  signal  processor  now  determines  when  to  perform  this  operation. 

Its  criterion  is  exactly  the  same,  i.e.,  reasonable  uniformity  of  echo  strength 
across  the  .rray.  Furthermore,  it  can  apply  an  additional  and  highly  sensitive 
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measure:  The  correlation  time,  which  is  the  reciprocal  of  the  echo  modulation 
bandwidth,  varies  inversely  with  target  size.  When  the  highlight  disappears 
the  "target"  in  (8)  becomes  the  size  of  the  aircraft.  If  the  one  meter  flat 
plate  reflector  were  on  a  30  m  aircraft  the  modulation  rate  would  jump  by  a 
factor  of  30  following  the  1/2  s  phase  synchronizing  period.  Similarly,  the 
echo  modulation  bandwidth  would  drop  by  this  factor  when  a  stable  echo  appears, 
thereby  permitting  the  signal  processor  to  readily  determine  when  to  phase 
synchronize. 

Observation  of  the  drop  in  echo  modulation  bandwidth  can  be  made  in  the 
low-frequency  control  branch  of  the  phaselock  loop,  if  such  a  circuit  is  used 
for  adaptive  beamf  rming,  or  in  the  data  store  of  a  real-time  digital  processor. 
In  the  latter  case  the  data  may  be  organized  in  the  same  format  shown  in  the 
lower  part  of  Figure  2  except  that  the  horizontal  coordinate  is  no  longer  range, 
but  time.  All  the  data  come  from  the  same  range  bin  and  the  columns  are  spaced 
by  one  interpulse  period.  Successive,  columns  are  successive  range-gated  samples 
of  the  echoes  at  each  antenna  element.  The  processor  observes  the  data  band¬ 
width  and  the  amplitude  stability  to  determine  when  to  phase  synchronize 
(Step  3).  During  the  dwell  time  of  a  highlight  the  signal  processor  ceases  the 
imaging  process  and  phase  conjugates  instead  (Step  A).  Next  it  observes  when 
the  modulation  bandwidth  returns  to  normal,  at  which  time  it  begins  again  the 
imaging  process  (Steps  7  and  8).  Following  the  high-resolution  angle  scan  it 
awaits  the  next  synchronizing  interval  to  repeat  the  process.  The  cycle  time 
is  the  order  of  one  second.  During  several  such  seconds  the  aircraft  has  moved 
and  rotated  insignificantly  insofar  as  the  observer  is  concerned,  but  has 
altered  its  orientation  sufficiently  to  cause  successive  images  to  be  statis¬ 
tically  independent.  Hence,  the  logical  next  step  is  to  form  "multiple  expo¬ 
sures,"  i.e.,  sums  of  intensities  of  several  successive  images.  Doing  so  builds 
up  the  information  content  in  the  image  and  fills  in  the  speckles  or  glint 
resulting  from  the  narrowband  (near  monochromatic)  radiation.  In  addition,  it 
reduces  the  sidelobe  peaks  of  the  random  array  by  several  dB  [1], 

An  experiment  has  been  conducted  to  measure  highlight  dwell  times  and  the 
percentages  oc  time  that  aircraft  targets  are  in  the  so-called  "synchronizing" 
and  "imaging"  modes  [10].  The  basis  for  the  measurement  is  the  expected  change 
in  the  probability  density  function  (pdf)  of  the  echo  amplitude  when  a  high¬ 
light  occurs.  Without  it  the  echo  is  the  vector  sum  of  a  large  number  of  small, 
randomly  phased  echoes.  The  sum  is  a  two-dimensional,  zero  mean,  random  Gaussian 
variate,  the  amplitude  of  which  is  Rayleigh  distributed,  i.e.,  if  A  is  the  echo 

2 

amplitude  and  2o  is  its  mean  square  value,  the  pdf  of  A  is 

f (A)  »  (A/o")exp(-A2/2c2) . 

When  a  strong  scatterer  of  strength  A  is  added,  the  pdf  changes  to  the  Ricean 
distribution  0 

f  (A)  -  (A/.72)exp(-(A2+A^)/2o2]Io(AAo,a2), 

where  I  is  the  modified  Bessel  function  of  the  first  kind  and  zeroth  order 
(which  reduces  to  the  Rayleigh  pdf  when  Aq  ■■  r)  [11]. 

The  experiment  was  designed  to  gather  simple  statistics  on  the  short-time 
pdfs  of  radar  echoes  from  aircraft.  The  radar  was  an  L-band  AN/TPS-1D.  It  is 
located  at  the  Valley  Forge  Research  Center  of  the  Moore  School  of  Electrical 
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Engineering,  University  of  Pennsylvania,  at  Valley  Forge,  Pennsylvania,  40  km 
west  of  Philadelphia.  Commercial  aircraft  flying  in  the  Philadelphia-New  York- 
Scrantou  triangle  were  observed  and  tracked.  Distances  varied  from  15  to  100 
km.  Echo  sequences  of  varying  numbers  of  pulses  were  sorted  by  amplitude  in 
real  time  in  a  histogram  generator  and  the  nature  of  the  pdf  was  judged  to  be 
Rayleigh,  Ricean  or  nondescript. 

The  results  are  shown  in  the  tables  below  for  broadside,  nose  and  tail 
aspect.  Sample  histograms  are  shown  in  Figure  10.  As  expected,  Rayleigh  and 
Ricean  pdfs  did  occur,  the  former  with  about  three  times  the  frequency  of  the 
latter  (Table  2).  About  20%  of  the  sample  runs  showed  no  preference  for  either, 
exhibiting  instead  tendencies  toward  uniformity  or  bimodality,  for  example. 

The  average  correlation  time  ir.  each  mode  also  was  measured  (Table  3).  Note 
that  the  1/2  s  dwell  time  estimated  earlier  for  the  sync  mode  is  within  the 
range  shown.  The  measurements  confirmed  the  expectation  that  aircraft  will 
alternately  permit  array  synchronizing  and  imaging. 


relative 

FREQUENCT 

OF 

OCCURRENCE 


RICEAN  PDF 
(A0*  So-) 


Fig.  10  —  Histograms  of  aircraft  echo  sequences 


Table  9  —  Relative  Frequencies  of  Occurrence  of  the  Imaging  and 
the  Phase  Synch)  onizing  Modes 


View 

Number 
of  Aircraft 

Imaging  Mode 
(Rayleigh  pdf) 

Synchronizing 

Mode 

(Ricean  pdf) 

Neither 

Broadside 

3 

.594 

.286 

.120 

Nose 

2 

.501 

.363 

.136 

Tail 

4 

.668 

.048 

.284 

1 


9 


.606 


.193 


.201 


Table  3  —  Approximate  Correlation  Times  for  Aircraft  Echos 


View 

Number  of 

Aircraft  Observed 

Average 

Correlation 

Time (seconds) 

Average 

Deviation 

(seconds) 

Broadside 

16 

1.56 

0.68 

Nose 

7 

0.30 

0.04 

Tail 

4 

0.17 

0.11 

SUMMARY 

A  radar  antenna  array  distorted  beyond  the  normal  tolerance  of  about 
one-tenth  wavelength  can  be  made  to  function  as  a  diffraction-limited  imaging 
aperture  by  adaptively  phase-compensating  for  the  distortion.  An  external 
point  source  of  radiation,  called  the  phase  synchronizer  or  adaptive  beamformer, 
illuminates  the  array.  The  source  may  be  an  active  beacon  or  a  passive  reflec¬ 
tor  echoing  the  radiation  from  a  radar  transmitter.  The  phase  of  the  radiation 
field  is  measured  at  each  array  element.  Phase  shifts  are  added  to  each  ele¬ 
ment  channel  to  eliminate  the  phase  differences.  An  array  focused  at  the  beam¬ 
forming  source  results.  Open  loop  scanning  in  range,  and  angle  follow  adaptive 
focusing  of  the  array. 

An  algorithm  suitable  for  digital  signal  processing  is  given.  It  describes 
(1)  a  search  procedure  for  locating  the  target  most  favorable  for  adaptive 
beamforming,  (2)  the  adaptive  beamforming  process  and  (3)  range  and  angle 
scanning  of  the  focused  beam. 

Experimental  evidence  of  the  validity  of  the  technique  is  given,  based 
upon  experiments  with  a  27  m,  X-band,  random,  sparse  array. 

A  modification  to  the  basic  algorithm  is  described  to  accomodate  an 
isolated  target  such  as  an  aircraft  or  a  ship.  A  radar  experiment  with  aircraft 
of  opportunity  disclosed  that  the  fluctuating  properties  of  airborne  target 
echoes  satisfy  the  requirements  of  the  search  procedure  for  a  suitable  adaptive 
beamforming  source. 
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